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A Computer Subroutine for Stress 
^_nalysis of Sotati.rg Disks - II. 

by John £. Brock 

Based upon theory developed by the '/vriter, R. E. Brown developed 
a successful computer program for analysis of radial and circumfer- 
ential stresses in rotating axisymmetric disks of variable thickness 
having an axisymmetrical thermal strain field. The '^ricer revised 
3rc’.'/n's program so as to Lnvoke a group of ancillary subroutines 
v/hich have been found useful in another application. In doing so, 
hov;ever, much unnecessary and copifusirig normalisation '.vas irxrcducec. . 

In particular, one of the normalisations would cause she analysis so 
fail in the quite ccmm.on case of a disk with no radial loading at its 
outer boundary. All this material appears as Reference 1, hereof. 

Referees evaluating a paper based upon Reference 1, called 
attention to these faults so t.hat the program has been revrritten. 

A listing of the main subroutine, RODISK, as revised, as 'well as 
listings of the ancillary subroutines may be found in Appendix A 
hereof. The reader will nose that other changes have also been made 
resulting in somewhat more flexibility of application. Employment of 
the revised program is described in the textual material which appears 
at the begirnir.g of the listing. 

Appendix B contains a revision of the second illustrative examole 
problem of Reference 1. This problem ',^/as solved for various values of 
= N-1, the number of equal subdivisions i.nto which the annular radius 
b-a is divided for purposes of numierical analysis by RODISK. Also, a 



number of different values of K?(3) were used. If KP(3) > 0, its 
value is the number of iterations which will be performed by RODISK. 
If K?(3) < 0, iteration will continue until tbiree successive values 
of the unknown parameter 3 decermined in the course of the analysis, 
satisfy the relation 

1 El - 82 1 + I 32-33 1 + IBi-Bi! 

< 10 '^^^^ 

iBiI + IB 2 ! + |B3| 



Vie also determined execution time by use of the library subroutine 
E<CL0K, executir.g under C?-cms on the IBM 36O/67 at the VI. R. Church 
Computer Center at the Maval Postgraduate School. 

We found that execution time per iteration is 

= 1.2 M + 5 (milliseconds) 



iter 



for any problemi. 



Accuracy vas evaluated by dealir.g with problems raving available 
aral:.n:ic solutions. It was found that the principal limitation on 
accuracy is determined by the choice of subdivions, the Integer M = 

N-1, so that there is a certain irJierent error regardless of how 
many iterations are made. This error depends on M, of course, and upon' 
the details of the problem. 'The error is greatest near the Inner ' 
radius of an annular disk, and is large if the ratio a/b is small. 
Fortuitously, the error may be smaller for an early iteration than 
for a somewhat later iteration but this is not practically useful 
information. For the problemi of Appendix 3 hereof, with a/b = .165, 
we find the results given in Table 1, (see next page). 

Thus, for example, with = 20, there is an irnerent error of 
about 1% and the results are not sigraf leant ly Improved by iteratir^ 



_ 0 _ 
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approx . 
limiting 
% error 


approx, 
iters, 
reo 'd . 


total 
time, 
secs . 


5 


16 


5 


.055 


10 


5 


7 


.12 


20 


1 


11 


.32 


40 


.1 


17 


.90 


100 


.01 


25 


3.1 


Table 1. 
execution 


Perce.nt error, required iterations, ai 
time for oroblem of Aopendix 3. 



r.ore than eleven tLT.es. V/ith eleven iterations, the solution 
is returned from RODISK in 0.32 seconds. 

The significant co.nclusion is that the execution is so fast 
that one rray as well take M = 100 (corresponding to N = 101, the 
rraximum available under present dimensioning) and iterate many 
more times than is strictly necessary. Takiag M = 101 and 
KP(3) = S gave execution in 3-7 seconds with 31 iteratiora and 
with an accuracy of 0.00^% (In the problan at hand, a^(a) 
specified as zero and the program gets -l.l^E-ll so that the error 
here is "infinite'.' Our evaluation of 0.00^% is for the first 
position rather than for the zeroeth.) 

This concludes the text proper of the present report. h^''/ever, 
we take advantage of this opportunity to correct errors in Refere.nce 
1 , viz . ; 

(1) Page 3, equation 12 should read 

.m = i/(n^-^vn+^) = ±/[ (n-2)^+^(l-v)n] 
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(2) Page 6, line 2. In place of T read aT. 

(3) Page 6, equation 33. Lower Unit of integration should 
be a rather than 0. 

(4) Page 7, line followi.ng equation 4 q, Reference should be 
to equation 37 rasher than eq'uation 33 - 
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Appendix A 



Listing of 
Subroutine 



RODISK 

and ancillary 
subroutines 



C SUBROUTINE POOISK. JOHN E. EROCK, 1 NAY 1978, PEVISEO 1 AUGUST 1978. fCDCCClC 
C This IS A SUBROUTINE ECR 0B7ERNIMING RAOIAc ANO CIRCUMFERENTIAL S7KESSRGCCOOSO 
C I N AN AXI SYMMETRIC THIN ELASTIC DISK HAVING AN AXISYMMETRIC THERMAL RQ0CCC2C 
C STRAIN FIELD AND ROTATING AT ANGULAR VELOCITY ONEGA ( RADI AN S/3 ECGNC ) RC0CCC40 
C ABOUT THE AXIS OF SYMMETRY. TWO TYPES OF FFOBLEM MAY BE TREATED: R0C00050 
C TYPE 1: ANNULAR DISK GF INSIDE RADIUS ARAD AND OUTSIDE RADIUS RCCCCCtC 
C BRAD. THE RADIAL STRESS IS SPA AT THE INNER RADIUS AND RCCCCC7C 
C SR3 AT THE CUTER RADIUS. THc INSIDE RADIUS MUST BE RCDCCCEC 
C GREATER THAN ZERO ROCOOCSO 
C T-YFE 2: SCLIC DISK HAVING RADIAL STRESS SRB ATOt/TSIOE RADIUS BRAD. PCOOOlOO 
C THE USER MUST PROVIDE A t'AlN PROGRAM VHICH CALLS SLBROUTINc RODISK RCDCCllC 
C AFTER IT HAS SLPPLIEC THE fCLLCWING INFCRmaT ICN, FCDCC120 
C m N, INTEGER. ( iS- 1 ) IS THE NUMBER Or EQUAL SUBDIVISIONS INTO Wh I CHRCD CO 1 20 
C THE ANNULAR RADIUS (oFAO MINUS ARAD) IS 0 W ICcC FCP CCMPU- ROCCGl-^O 
C TATIQNAL FURPCSSS. THE PRESENT OIHENSICNING CAfl ACCCMHQDATE N RCDCClfC 
C NOT GREATSR than 101. , ROOOOIEO 
C (2) BRAD RCDCC170 
C (3) ARAD (NOT NECESSARY FOR PROBLEMS OF TYPE 2.) RCCCCIEO 
C (.4) SRB POOCC19C 
C (5) SRA (NOT NECESSARY FOR PROBLEMS OF TYPE 2.) PCDCC2C0 
C (6) PGIS, POISSON'S RATIO FODC02 10 
C Ul KPlD^rl.E. INTEocP TO DENOTE PPCBLEM OF TYPE 1,2. RC0CC22C 
C (81 K?(2). INTEGER TO PROVIDE FOR SKIPPING I-rilLE PRINTING FCCCC22C 
C OUTPUT. FCR example, IF N=101 AND <F(2)=5, ONLY EVERY PCDCC24C 
C FIFTH SET OF VALUED isl LL BE PRINTED; 1ST,6TH,..., 96TH, RCC0O2fO 
C AND lOlST. FG0CC26O 
C (5) KP(3), INTEGER SPECIFYING NUMBER OF ITERATIORE TO 8c PGCCC2IC 
C PERFORMED. USUALLY NP(3) = 10 IS SUFFICIENT FOR ENGIN- POCCC280 
C EERI.N& ACCURACY . ALTERNATELY, IF KP(3J IS A NEGATIVE RCCC029O 
C INTEGER, ITERATICN WILL CONTINUE UNTIL THREE SUCCESSIVE PCDCC2CC 
C VALUES OF A PARAMETER E, DETERMINED INTERNALLY, ARE RCCCC2IC 
C SUFFICIENTLY CLCSE AS COMPARED TO AN ERSILCN EQUAL TO PCCCC22C 
C TEN RAISED TO THE KP(2) PQwER. RCCCC22C 
C (r-0) KP(4). IF KP(4)-0, GNLY FINAL ANSWERS WILL EE PRINTED, RCDCC2RC 
C IF KP(A)=1, A SEQUENCE OF ITERANT VALUES QF E WILL BE PCCCC25C 
C PRINTED TO INDICATE DEGREE CF CONVERGENCE. IF KP(R)>1, PCCCC26C 
C THERE WILL BE NO PRINTING AT ALL WITHIN RGOISK, BUT UPON RCCCC27C 
C RETURN KP(5) WILL CGNTAIN THE NUMEER OF ITERATIONS VHICH PCCCCEEC 
C WERE PERFGRMEO SO THAT KP ( 5 ) MUST BE RESET E.EFORE RODISK PCOCCSSC 
C IS CALLED AGAIN. RC0C04C0 
C (11) KP(5). KP(5I=0 CALSES MILNc CUBIC 'SPLINE IMEGRATICN PCCCC410 
C TO BE USED. OTHERWISE TRAPEZOIDAL INTSGRA7ICN IS USED, RC0CC42C 
C (12) VECTOR X(1,J), J=1,2,...,N, CONTAINING VALUES OF DISK FCDCC42C 
C THICKNESS AT EQUALLY SPACED RACII FRCM INSIDE TO OUTSIDE. PC0CC4RC 
C (13) VECTOR X(2,J) CONTAINS VALUES CF GAMMA TIMES ONEGA RCCC04SC 
C SQUARED WhIrE GAMMA IS (MASS) DENSITY Q= The MATERIAL. PCOCCAAC 
C FOR MOST FFuBLEMS GAMMA DOES NOT VARY WITH RADIUS AND RCCCC47C 
C ALL ELEMENTS CF THE VECTOR WILL BE THE SAME- RCDCCREC 
C (lA) VECTOR X(2.J) CONTAINS VALUES CF ( EE ) ( A L P HA ) ( T £ £ ) WHERE PCCCC49G 
C EE IS YOUNG'S MODULUS, ALPHA IS THE COEFFICIENT OF LINEAR RCCCC5CC 
C THERMAL EXPANSION, ANO tE£ IS THE T £M P B R A 7U RE CHANGE. PCDCCSIC 
C THE MAIN PROGRAM MUST CONTAIN THE STATEMENTS: FCCCC520 
C IMPLICIT REALMS (A-H,C-Z) • PCGCC52C 
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PEAL^S X ( 20t 101) 

INTEGER K?(5) 

COMMON 

CCMMCN /CNE/^=lAO.eR^iOtSRA.SPStPOIS 

FOLUG’^I.NG SUaRGi-TIiNc ROCISK THERE ARE SEVERAL ANCILLARY 
SU8RCUTINSS ^HlCr PERFORM vARiCUG CPEPATICNS CN THc VECTORS 
X(ItJ). THE PURPOSE OF EKH IS 06V ICU5 FRCf- THE LISTING. 

THEY MAY 8E EMOlQYED IN THE USE R * S MAIi\ PRCGRAM. SbERCUTINE 
OUPVf which CUPLICaTES a VcCTCRt AND SU8RCUHNE PRiVt «FICH 
PR.IMS A VECTCRt ARE NCT CALLED cV POOISR cLT MAY 6c USEFLL 
IN THE USER*S MAIN PROGRAM- 
SL8ROUTUE RCOISK 
I f'FLI Cl T RE AL=»6 ( A-HtC-Z ) 

REALMS X(20, ICl) 

INTEGER KP(5) 

CCMMGN XtNtKP 

CCMMGN /CNE/ARAC,SRAO,SRAt SRSt. PQI 5 

onE=i .o-»o 
ZERO=0-C^0 
61=1.0^6 
82^1.0+10 
63= 1 . 0+12 

I F (KP ( 3 ) •LToO ) EFS= ( 1 ) ^^(KP ( 3 ) J 

I fiKPi 1) .EO. 2) ARAC=ZSRC 
EMA= 6RA0-ARA0 
ENM=N- 1 

I f (jCP ( A ) .LE .1 ) V^RITE(6t2J KP(I) 

2 FCRMATt// t IDX ,*RGuISK PROBLEM OF TYPE 
CO 5 1=1 tN 

E1M= I- 1 
Y=EIM/ENM 

X (4,n =ARA04(8RAD-ARAC)^Y 
X (5tl ) =Y 
5 X (6tn ^Y 
I TER= 1 

IF( KP( 1).60.2J GO TO ICC 
: THE PROBLEM IS CFTYPc 1: ANNULAR DISK 



Cl= ( 2 -C^C+PQ I SJ 3R AD-ARAO J 
CALL INTV( 1 ,7 tBMA ) 

C2=X(7,N J 

C5=X( 3 , ] )-X(3 tN ) + < CNE- P 01 S)=5« (SR A-SRE) 
CALL MLL V<1 t2 r3 J 
CALL MULV(4,6,S) 

CALL INTV(9tlO,aMA) 

C6 = X(10tN) + X(itN) ^SR8-Xa» D^SRA 
20 CALL I NTV ( 6 rll , SMAi 

C3=6RA0+ (CNt+PG I S) ^X( 1 ItN) 

CAL L MU LV( 1, 6, 12) 

CALL I NTV( 12t lit 3MA) 

C^ = X( 13 tN) 

0=C1^C4-C2^ C3 
A=(C5^C4-C6^C3) /O 
8= (Cl ’>C6-C2^C5) /O 
30 CONTINUE 

IF (KP(4).EC.i) WRITE(^,T) ITERtAtS 
7 FCRMAT(5>t I10tlF2E20-5) 

CALL MULS(7t lAtA) 

CALL MULS (13t 15 t c ) 

CALL AOCV( IZt 15, 151 
CALL SU6V( 15t ID t 16) 

S=SR8^X ( 1, N )-X( 16 tN ) 

CALL AODS(UtlctS) 

CALL 0 IVV( 16t 1, 16) 

ZA=X (3 t 1) + A^'ARA0+ (CN E'-POIS ) ^X (16 , 1) 
CALL MULS(lltl7tB ) 

CALL sues (4 , 13 tARAC ) 

CALL MULSi 18tl8 ,A ) 

CALL ADOV( 17,13,1 1) 

S=-( GNE^PCIS ) 

CALL MUIS( 17, 17,3 ) 

CALL ADCS tl7 t 17 tZA) 

S=CNE-PC IS 

CALL MULS (16,13,5) 

CALL SU8V( 17,3, 17) 



RCCCC5A0 
ROOCC55C 
RC000560 
P00CC5 70 
RCDCC5 £C 
PC0CC5SC 
RCDCC6CC 
R0CG06 10 
PCDC06 20 
ROCGCdEC 
R0CC0640 
RaCCC65C 
RO CCC660 
RCC CC670 
ROCCG630 
RCCC06S0 
RCCCC7CC 
RCDCC71C 
RGCC0720 
RC0C07 30 
RGCCC7AC 
ROCC0750 
RCCC076G 
RO CCC77C 
RC0C0780 
RaOCC7SC 
RG0CC3C0 
RGCCCa 1C 
RCDCCS2C 
RCCC08ED 
PCDC0840 
ROCCC65C 
ROCCC66C 
RCCCC67C 
ROCCC88C 
PCCCC8FC 
RCCCCSCO 
RCCCC910 
RCCCC92G 
RCDCC93C 
RGDC09^0 
PC0C0950 
RCC CCS 6C 
R00CC970 
RCOCCSaC 
RO CCC990 
RCO CIOGO 
RGC Cl Cl C 
ROCCICZO 
RGCC1C3C 
ROOCICZC 
RCCC1C50 
PODC1C60 
RCDCIC TC 
ROCC 1C8Q 
RCDCIC^C 
ROC Cl ICO 
RCOClllO 
RCCC112C 
RCDC113C 
R0CC114C 
RO cell 5Q 
RCCC1160 
ROD01170 
RGCCll 80 
ROCC 1190 
RCD C12C0 
RCCG1210 
PC0C1220 
RC0C123C 
RCOC 12^0 
RCCC125 C 
ROCC1260 
RCC C1270 
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C4LL suev( 17 , 13 ♦!? ) 
nER = ITEP + l 

I MKP( 3 ) -LI.C ) GJ TG 215 
IFi ITER-GT.XP (3) ) GJ TG 200 
40 CCNT INUE 

C;iLL 0UPV(4tl^) 

I F (KP ( 1 ) ,50 .2) X( 14, 1 J = C^e 
CALL 01 VV ( 17, 1 ^,1 t) 

C^LL SUBS ( i£, IS, A) 

CALL 01 ^S I13,c,aj 
1F(KP( 1 ) -EQ.2 J GO TO 15C 
GC TO 20 
100 A=2SRC 

C=(ONE-POIS )/X( 1,1) 

CALL MULV(1,2,7 ) 

CALL MULV( 4,7,8) 

CALL INTV (8 , 10 , 8MA ) 

SUM = X ( 10 ,N ) +SR6^( X( 1 , N )-X ( 1, 1) ) 

SLJ^=C* SUM+ X ( 3,1 )-Xl3 , N) 

150 CALL I NTV(6,11,3MA) 

0EN = 6RAC+ {0^6 + P01S)^xal,N ) 

CALL MULV(1,6, 11) 

CALL 1NTV( 11, 13, 0HA ) 

D£N=OEN+C^X (13 ,N) 

8=SUM/0eN 

CALL INTVd ,7 ,8.‘^A > 

CALL INTV(6,11,8MA) 

GC TO 30 
215 B3=62 
B2=8 1 
81=8 

2UM=0AES( ei-62) +0 ASS ( c2- E3 )'*■ CABS ( 3 3-3 1 ) 

0IV=0A3S <61 )+0 A eS( 82) +CA3S( 83 ) 

CPIT=2UM/CIV 

I F(CRI T. LT-EPS) GC TO 200 
GO TO 40 

2C0 CALL ACCV(17 ,16 ,19; 

IF(KP(2) .LT-G.ANO .KP(4 J,£Q*0) '«R 113(6,201) ITER, EPS 

201 FORMAT (//, 20X ,18 , • ITERATICNS RECUiPcO \sllH EPSILCN = •,1PE8.1) 
IF(KP(4 ),L£*1 ) i.RIT£(6 ,2C4) 

204 FORMAK//) 

IF(KP(4 ) ,LE.l ) WRITS(6,2C5) 

205 format ( 23X , ’RAOIUS’ , 10>, *THI CKNES5 •, 5X, • GAMMA OMEGA SQ«,7X, 

I‘EE ALPHA TEE • ,7X ,♦ SI GMA RAD I AL * t 6X , • SI GM A CIRCUMFM 

NSKIP=KP (2 ) 

00 210 I=1,N,NSKIP 
J = I/NSKI F 

1 F(KP (4 ) •LE. 1) .VRIT E (6, 211 )J , X( 4, 1 ) ,X( 1,1) ,X( 2,1 ) , 

1X< I ) , X( 1 6,1 ) , X ( 19, I ) 

210 CCNT INUE 

I F(KP (4 ) ,G 1 . 1 ) KP (5 ) = ITER 

211 FCRMAT( no , iPoE 19. 5) 

RETURN 

END 

C THIS IS THE START OF THE ANCILLARIES 
SUBROUTINE ADDV ( N 1 , N2 , N3 ) 

REALMS X (20,101) ,S 
INTEGER KP(5) 

CCMMON X,N,KP 
DO 1 1=1, N 

1 X (N3, I J=X(N1, i)+X(N2, I ) 

RETURN 

ENC 

SUBROUTINE SU8V ( N 1 ,N 2 ,N : ) 

REALMS >(20, 101), S 
INTEGER KP(5) 

COMMON X,N,KP 
00 1 I- 1 ,N 

1 X (N3, I )=X(N1,I)-X (N2, I) 

RETURN 
6 NO 

SUSROUTINE MULViNl, N2, N3 ) 

REAL’i'8 >(20,101) ,S 
INTEGER KP(5J 



POCC 1280 
RCDC12<C 
R0CC13CG 
PCDC13 10 
RODC132C 
R0CC1330 
RGCC134C 
RO CC13 5 0 
RCC C1560 
ROCC1370 
PCOC1350 
RQCC 1390 
PCDC14CC 
R0D01410 
PCD01420 
RCDC143C 
FCDC 1440 
RCDC145C 
ROCC146C 
PCC C147C 
ROC014 K 
R0001490 
RCCC15 CC 
R0DC151C 
RCCC1520 
FOOC 1530 
RCDC154C 
RQCC 15 50 
RCOC156C 
ROCO 1570 
PC0C1530 
ROCCl 5 SC 
R0D016C0 
R0CC161C 
ROOC1620 
PCC C163C 
RODC1640 
RGCC165a 
R0DC166G 
RCCC16 7C 
RCC016E0 
PCDC1690 
RCCC17CG 
F00017 10 
RCCC172C 
ROCC173C 
RCC (174C 
ROCCl 750 
POCC176C 
PCrCITTG 
RCDC17EC 
RCCCl 790 
PCOC18CO 
RCDCIEIC 
ROCC1S20 
RCDCl 630 
ROCOIS^O 
PCDC1850 
RCOCl 66C 
RO0C16.70 
RCCCISEO 
ROOC1S90 
RCDC19C0 
PC0C19 10 
RCDC 192C 
FCCC 19 30 
RCDCi S4C 
ROCC1950 
PCC C1960 
R0DC1S7C 
ROOC1960 
ROC C199C 
ROOC2CCO 
RCD02010 
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CCMMON XtNtKP 
CC 1 1 = 1, 

1 XUO .U (N 1, 1 (N2 , I ) 

RETURN 

END 

SUBROUTINE C I VV ( N 1 , N2 ,N 2 ) 

RE;iL^8 X(20,101),3 
INTEGER KP(5) 

COMMON X,N,XP 
CO 1 1=1, N 

1 X (N3 , I J =X( M T I i /X (N2 , I ) 

RETURN 

END 

SIBROUTIN5 /DCS(Ni ,N2 ,S) 

REALMS X(20,101)tS 
INTEGER KP ( 5) 

COMMON XtNtKP 
DO 1 1 = 1, N 
1 X( N2 ,I ) => (N1 tl ) ^5 
RETURN 
END 

SUBROUTINE SUB S ( N 1 ,N2 , S ) 

REiiL>5^8 X(20,101),S 
INTEGER KP(5) 

COMMON XtN,KP 
DO 1 1=1, N 
1 X {N2 , I } =X{ M t I }-S 
RETURN 
END 

SUBROUTINE MULS (N1 , N2 , S ) 

REAL=^8 X(20,101},5 
INTEGER KP ( 5) 

COMMON >tNtKP 
DO 1 I=1,N 
1 X(N2tI ) =X(.Nl,Ii^S 
RETURN 
END 

SUBROUTINE C I VS ( N 1 , N2 ,S ) 

RE;iL^8 X(20ti01),S 
INTEGER <<P (5 ) 

COMMON XtNtKP 
DO 1 I=ltN 
1 X (N2 , I)=X( N1 tl ) /S 
RETURN 
END 

SUBROUTINE PRIV(N1,I,JJ 
RUL^e X(20,lJi)tS 
INTEGER KF< 5 ) 

CCMMCN >tNtKP 

IF SECOND /5RGUMENT EC'JALS Ot GO DIRECTLY TO RETURN 
IF SECCNO ArGLMENT ECJALS 1, PRINT The VECTOR. 

IF SECOND ARGUMENT EC’JALE 2, PRINT THE IDENTITY AND THE VECTOR. 

IF SECOND ;iRGLM2NT ECJAL3 3t PRINT THE VECICR NUMBER ANC THE VECTOR. 

IF SECOND ARGUMENT ECUAUS ^t PRL'U NUMBERt IDENTITY, ^NC VECTOR. ' 

IF SECOND ;5RGUMENT EQUALS f, PRINT IDENTITY ONLY . 

IF( I .EQ .0) GO TQ iO 
I F( I .£Q .1 } GO TO 1 
IF( I .EG .2 ) GO TO 2 
I f( I .Eg .31 GC TO 3 
IF(I.EC.A) GO TO A 
IF( I .EQ .5) GC TG 5 

1 DC 8 K=1,N 

8 'rtP ITE (6 tS ) J<,X( M ,K I 

9 FCRMAT (30X, I5tiP62C.5 i 
10 RETURN 

2 WR1TE{6,21) J 

21 FCRMAT(//,30X,« VECTOR «]TH IDENTITY *,I5t* FOLLOWS ^ M 
GC TO 1 

3 WFITE (6,31) N1 

21 FCRMAK //t30X, « VECTOR NUMBER ‘,15,* FCLLCWS^M 
GC TO 1 

A HR ITE(6,41 } Mt J 

41 FORMAK // ,20X ,» VECTOR «ITH ICENTITY * , 15 , ‘ FXLLCHS:*} 

GC TC 1 



ROCC2C20 
R0CC2C3C 
ROCC2C4C 
RCCC205Q 
RODC2CEO 
RCDC2C TC 
P0DC2C80 
RCO C2CSC 
RCCC21C0 
RCDC21 10 
RCDC212C 
R0CC2120 
RCCC214C 
ROCC2IBO 
RCCC21EC 
R000217D 
RQDC21E0 
RGCC21S0 
PODC22CC 
RCCC2210 
FODC 2 2 2J 
RCDC223C 
POCC22AO 
PCDC225C 
RCCC2260 
FC0C2270 
RCCC22 U 
ROCC 2290 
RCCC23CC 



ROCC221C 



PCCC232C 
R000233D 
R0DC234C 
RCCC23 50 
FODC23EC 
RCCC23 TC 
ROD023cO 
RCDC23SC 
RCCC24C0 
RCDC24 iC 
R0CC2420 
FCDC24 30 
RCDC24^C 
ROCC245C 
RCCC24EC 
ROCC247C 
FCCC24EC 
ROC024SO 
RODC25CO 
FGCC2510 
RO D C 2 5 2 C 
FCCC252G 
FCDC25nO 
RCDC255C 
FCDC2580 
PCDC25TC 
RCCO 25 EO 
PCDC2590 
RCDC26CC 
FCD026 10 
RCCC262C 
ROCC262C 
9CCC2t^.C 
ROCC2650 
R0DC26E0 
RGCC2670 
POCC26EG 
PCCC2690 
F00027C0 
RCDC2T1C 
PCCC272C 
PCDC2T3C 
RCCO 27 <0 
FC0C2750 



I 



8 



i 



c 

c 



5 WPITE (6t51) J 

51 f CRMAK //,30Xt* VECTOR ^ITH lOcNT ITY ‘tlS,* HAS EEEN GENERATEC,*) 
GO TO 10 
END 

SUBROUTINE CUPyiNl»iMZ) 

REAL^J'B X(20,101},S 
INTEGER KP(5) 

CCMMON X,N,KP 
DO 1 I=ltN 
1 X (M2»I )=X(.N1, n 
RETURN 
END 

SU6RGUT INE I*NTV(N1,N2,S) 

KP{5}=C CAUSES MILNE INTEGR-^^T IGK TO BE LSED^ 

OTHERWISE TRAPEZOJCAL INTEGRATION IS LSEO, 

REAL* 8 X{ 20, ICl Sf AOC ,N INO ,NTN(j ♦ F IVQ ,ThT C »EN , R» F 
INTEGER KP{5) 

CCMMON X^NtKP 

I F(KP( 5) .NE ,0 ) GO TO 10 

EN= N-1 

EN=1 ,0+0/EN 

MNCh=EN^^,C+0/2. AO+I 

NTNO=eN*l,90+l/2. 

F I VO=EN*^ 5. Q + 0/2* AC + 1 
THTO=EN^1.30^'l/2.sC+i 
R^£N/2. AO+1 
X(N2, 1 ) = 0.0 + 0 

X (N2 , 2) =M-NC’J'X (Ni,l) + NTNC^X(NI,2)-FIV0*X(Nlti)+X(M,^)*F 

NM3=N-3 

CO 1 K^1,NM3 

KP1 = K + 1 

KP2=K+2 

KP3=K+3 

ACC-THTO’<«( X (N1 ,KP1 ) + X (N1 tKP2 n-R’^(X(NltKMX(Nl,KP3n 
1 X{N2,KF2)=X(N2,KPl)^‘Ano 

X (N2,N) =X( N 2, N-i ) ♦fIN5^X(Nl,N)+NTNC^X(Nl ,N-1)-FIV0^) (NltN-Z) 
1+XNl ,N-3)*R 
CALL MLLS(N2,N2,S) 

RETURN 

10 CCNTINUE 

X (N2, 1 )^0-C+0 
P=2^(N-1) 

DC 11 1^2, N 
J=I-1 

11 X (N2, I ) =X(N2 t J) {N1 t I }/F^=X(Nl, J J/F 
CALL MULS(N2,N2tS) 

RETURN 

END 



PCC C27CC 

PCCC271C 

PCCO 2760 

PCOC27SO 

RCDC2 ECC 

PCCC28 10 

RC0C2 E2C 

ROCC28SO 

PC002SA0 

RC0C2650 

RO0C2860 

PCCC267C 

ROOC26EC 

RC002890 

FOOC29CO 

RCDC2S1C 

POC C29 20 

RCD C2S3C 

RCCC 29 iQ 

PCDC2950 

RCC C2 <tC 

ROCC 2970 

RCC C2 SEC 

ROCC2990 

RCC C3CC0 

ROCC3010 

RCCC3C20 

P0CC3C2C 

RCDC3C AG 

RGCC3C50 

PCDC3C60 

RCCC3C iC 

ROCC3CcO 

ROCC3GSO 

PCDC3 ICO 

RCCC311C 

PCDC3 1 2C 

PCCC3 13C 

PCDC31 

PCCC3 ] 5C 

PCCC326C 

PCCC317C 

PCCC3 16C 

PC0C319C 

RCDC32CC 

PCCC321G 

RCCC322C 

FCDC323C 
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Appendix 3 



Sample Problem 

A disk rotating as 7200 rpm and composed of a metal caving a 

a specific weight of 0.283 pounds per cubic inch, is 0.35 inches 

inside diameter and 5.15 inches outside diam.eter. The radial stress 

at the inside radius is zero and that at the outside radius is 22,000 

psi. The thickness ’/aries with radius according; to the law 

—0 42 

t = 0.1^93 r"'^* ' (all dimensions in Inches) 
and the temperature change (from the zero stress condition) is given by 
T = 60 - 1.6 r^. 

Take E = 29,000,000 psi and a = 6.7 10”° /°? and determine radial 
stress (a^) and circumi’erential stress (a^ ) as f’onctions of r. 

This problem, illustrates most of the capabilities of RODISK, 
Because of the special nature of the thickness ’/ariation, i.e, a 
power relation, an analytic solution may be established so that 
the accuracy of the RCDISX solution may be evaliated. Results of 
such evaluations are given In Table 1 of the body of this report. 

There it may be seen t.hat accuracy far better than en.gineerin.g 
considerations require or justify may be obtained by taking, say, 

N = 101 and XP(3) = 25, so that in 3.1 seconds RODISK returns to 
the calling (i.e., input) program results with a maximum error of 
0.01 % or less. The tabulation which follows shows output with N 
= 101 and KP(3) = -6, resultin.g in 27 iterations and taking 3.3 
seconds. Accuracy is better than . 006 %. 
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the output shown here, Execution tiiae was 3.3 second; 
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